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REDUCTION  OF  BEAM  BREAKUP  GROWTH  BY  CAVITY 
CROSS-COUPLINGS  IN  RECIRCULATING  ACCELERATORS 

I.  Introduction 

1  2 

The  spiral  line  induction  accelerator  (SLIA)  ’  is  a  novel  compact 
accelerator  with  the  potential  of  accelerating  a  multi-kiloamp  electron 
beam  to  tens  of  MeV  or  beyond.  It  makes  use  of  a  stellarator  magnetic 
field  to  force  the  beam  to  recirculate  through  the  same  accelerating  unit, 
and  is  projected  to  share  many  of  the  advantages  of  linacs  and  cyclic 
accelerators.  A  proof-of-concept  experiment  (POCE)  is  currently  under 
construction.  It  is  designed  to  trap  a  35  ns,  10  kA  beam  and  accelerate 
it  to  8.5  MeV  through  two  recirculations.  A  schematic  drawing  of  the  POCE 
is  shown  in  Fig.  1.  It  consists  of  two  accelerating  units.  The  beam  path 
threading  through  these  units  consists  of  four  arms,  labeled  sequentially 
as  A,  B,  C,  D.  Each  of  the  four  arms  consists  of  five  accelerating  gaps; 
there  are  twenty  gaps  altogether.  The  gaps  in  Arm  A  may 
electromagnetically  be  coupled  to  the  gaps  in  Arm  C  and  those  in  Arm  B  to 
those  in  Arm  D. 

Since  beam  acceleration  occurs  only  when  the  beam  encounters 
accelerating  gaps,  attainment  of  high  energies  necessarily  requires 

4-10 

sequential  beam  passage  through  these  gaps.  Beam  breakup  (BBU)  growth 
is  then  a  primary  concern  of  the  viability  of  the  SLIA  configuration  to 
achieve  ultra-high  energies.  In  addition,  it  is  necessary  to  devise  a  new 
procedure  to  assess  the  BBU  growth  in  the  POCE  since  scalings  for  such  a 
configuration  are  unavailable. 

The  SLIA  is  a  hybrid  between  a  linac  and  a  cyclic  accelerator.  BBU 
growth  would  occur  in  the  linear  sections  of  the  accelerating  units  where 
the  gaps  lie.  However,  since  the  beam  recirculates  through  the  same 
accelerating  units,  albeit  via  different  beam  paths  (Fig.  1],  there  would 
be  inevitable  cross-coupling  of  the  gaps.  Thus,  the  first  gap  of  Arm  A 
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would  couple  CO  the  first  gap  of  Arm  C,  the  second  gap  of  Arm  A  would 
couple  to  the  second  gap  of  Arm  C,  etc.,  {Fig.  2].  In  fact,  there  has  been 
considerable  concern  that  such  cross-couplings  would  worsen  BBU  growth. 

This  concern  has  prompted  consideration  of  cavity  designs  that  would 
minimize  such  cross-couplings.^^  (We  will  assume  here  that  the  beam  pipe 
is  cut  off  to  the  propagation  of  the  cavity  modes  responsible  for  BBU,  so 
that  the  couplings  illustrated  in  Fig.  2  are  the  only  ones  possible.) 

Since  the  SLIA  is  markedly  different  from  a  linac,  an  attempt  was  made 
to  assess  its  BBU  growth  by  taking  the  opposite  view,  pretending  the  SLIA 

9 

to  be  a  strictly  cyclic  system.  That  analysis,  though  not  directly 

applicable  to  SLIA,  does  pinpoint  some  intrinsic  differences  between  BBU  in 

9  10 

a  linac  and  in  a  cyclic  accelerator.  ’  It  also  raised  several 

interesting  Issues  concerning  whether  BBU  growth  in  a  cyclic  system  is  best 

described  in  terms  of  transient  amplifications  or  in  terms  of  unstable 

eigenmodes,  as  the  properties  of  growth  depend  very  much  on  the  pulse 

length,  on  the  number  of  pulses,  on  the  pulse  separation,  on  the  quality 

factor  Q,  and  on  the  degree  of  feedback.  Many  of  these  issues  remain 
12 

unsettled. 

In  this  paper,  we  treat  BBU  in  a  SLIA  under  more  realistic  assumptions 
and  employ  parameters  that  cover  the  range  of  the  POCE.  The  crucial 
feature  of  this  work,  which  is  different  from  all  others,  is  the  explicit 
modeling  of  cross-coupling  of  cavities  in  the  different  arms  of  SLIA  (Figs. 
1  and  2].  Cavities  within  the  same  arm  are  shielded  from  one  Piiother 
electromagnetically ,  as  in  the  POCE.  Thus,  we  may  unwrap  the  SLIA  into  a 
linac  to  treat  BBU,  [Fig.  2],  except  that  cavity  number  1  is  coupled  to 
(and  only  to)  cavity  11,  cavity  number  2  to  number  12,  etc.  (The  twenty 
cavities  are  numbered  sequentially  in  the  order  of  beam  encounter). 
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The  moslr  unusual  result  that  ve  found  is  that  the  presence  of  cross¬ 
coupling  may  actually  reduce  BBU  growth.  Several  reasons  may  be  given  to 
explain  this  unexpected  phenomena:  the  relatively  short  pulse  length 
treated,  the  sharing  of  the  mode  energy  by  another  cavity  (and  therefore 
the  reduction  of  beam  deflection  by  an  individual  cavity),  and  the  phase 
mixing  associated  with  the  detuning  in  the  breakup  mode  frequencies  as  a 
result  of  finite  (though  small)  coupling.  This  reduction  of  growth  is 
observed  regardless  of  the  value  of  Q  associated  with  the  deflecting  dipole 
mode.  Based  on  this  study,  we  conclude  that  BBU  in  the  POCE  of  SLIA  is  not 
likely  to  cause  beam  disruption. 

The  following  assumptions  have  been  made  to  reach  the  above 

conclusions.  All  gaps  are  identical,  in  the  sense  that,  when  isolated, 

each  admits  a  deflecting  mode  with  the  same  traverse  shunt  impedance  Zj^  and 

the  same  quality  factor  Q.  Ve  assume  that  cross-coupling  would  not  alter 

these  values.  The  beam  transport  between  two  successive  gaps  (including 
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that  around  the  bend)  is  modeled  by  2  x  2  matrices  '  with  constant  matrix 
elements.  Ve  consider  only  a  continuous  coasting  beam,  and  assume  that  the 
beam's  pulse  length  is  less  than  1/2  of  the  total  transit  time  from  the 
first  cavity  to  the  last  cavity  so  that,  at  any  moment,  the  beam  passes 
only  one  cavity  in  the  cross-coupling  interaction  [Figs.  1,  21.  Ve  model 
an  individual  cavity  by  an  RLC  circuit  whose  natural  frequency  and  Q 
coincide  with  the  respective  values  of  the  deflecting  mode.  The  cross¬ 
coupling  between  the  cavities  is  modeled  by  a  mutual  inductance  linking  the 

equivalent  RLC  circuits.  It  is  done  in  much  the  same  way  as  cavity 

13 

coupling  is  handled  in  standard  microwave  literature.  The  formulation  is 
therefore  consistent  with  the  established  results  in  the  appropriate 
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limits:  the  BBU  growth  reduces  to  that  of  the  cumulative  BBU  in  a  linac^’^ 
when  the  cross-coupling  is  absent  and  the  standard  coupled  cavity  result  is 
recovered  when  the  be2un  is  absent. 

In  Section  II,  ve  describe  the  cross-coupling  of  the  deflecting  modes 
via  simple  RLC  circuits.  There,  we  introduce  the  dimensionless  cross¬ 
coupling  coefficient  K.  In  Section  III,  we  formulate  the  governing 
equations  for  BBU  evolution,  including  the  effects  of  cross-coupling, 
together  with  the  imposed  boundary  conditions  and  initial  conditions.  In 
Section  IV,  we  present  the  numerical  results  over  a  wide  range  of  K  and  Q. 
It  is  these  numerical  results  on  which  ve  base  our  conclusions,  stated  in 
the  Abstract.  In  the  last  section,  ve  discuss  further  the  nature  of  cross¬ 
coupling,  and  examine  the  possible  consequences  when  ve  relax  some  of  the 
simplifying  assumptions.  The  implications  are  explored.  Ve  also  raise  a 
number  of  issues  that  need  to  be  resolved  if  the  SLIA  is  to  be  scaled  to 
much  higher  energies  than  the  POCE. 


4 


II.  Cross-Ccapling  Between  Cavities 

When  the  gaps  are  isolated  from  one  another,  ve  assume  that  the 

deflecting  mode  may  be  modeled  by  an  equivalent  RLC  circuit  of  frequency 

and  quality  factor  Q.  Vhen  cross-coupling  between  two  gaps  (cavities)  is 

present,  these  individual  RLC  circuits  are  also  coupled.  The  coupling  of 

13 

cavities  has  been  customarily  modeled  with  a  mutual  inductance  M  e  kL 
which  links  the  RLC  circuits  representing  the  individual  cavities  [Fig.  3]. 
Here,  K  is  the  dimensionless  constant  which  measures  the  degree  of  (cross-) 
coupling.  Vhen  the  coupling  is  weak, 

K  «  1.  (1) 

Ve  assume  that  the  presence  of  cross-coupling  does  not  alter  the  values  of 

Q  and  of  the  deflecting  mode  in  an  isolated  gap.  Thus,  =  1//LC  and  Q 

=  to  RC. 
o 

To  account  for  the  cross-coupling  between,  say,  cavity  1  of  arm  A  and 
cavity  11  of  arm  C,  [Fig.  2],  we  first  denote  /^^^t)  to  be  the  deflecting 
mode  amplitude  of  the  i-th  cavity  at  time  t.  Vhen  the  weak  coupling 
condition  (1)  is  satisfied,  and  may  be  shown  to  be  governed  by 

the  following  differential  equations  when  the  beam  is  absent  [see  Appendix 


A]. 

(2) 

(3) 

where  L  denotes  the  operator 

L  a  d^/dt^  +  (Wjj/Q)  d/dt  +  </. 

(4) 
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Equations  (2)^  and  (3)  admit  normal  solutions  of  the  form  exp(ia)t).  When  K 

=  0  (zero  coupling),  both  and  oscillate  independently  with 

2 

exp(i(i)^t  -  (0^t/2Q)  dependence.  Vhen  K  «  1,  the  eigenfrequencies  are 
given  by 

“  =  Ho  •  <1±K)]<.„.  (5) 

Equation  (5)  shows  that  the  coupling  leads  to  a  slight  detune  of  the 
breakup  mode  frequency  and  that  there  are  two  independent  modes  of 
oscillations. 

The  coupling  between  cavity  number  2  and  number  12  is  also  described 
by  Eqs.  (2)  and  (3),  in  which  is  replaced  by  and  by 

A  similar  procedure  applies  to  coupling  between  cavity  number  3  and  cavity 
number  13,  and  so  on. 

We  have  assumed  that  there  is  no  phase  shift  in  the  coupling  between 
cavity  number  1  and  cavity  number  11.  If  the  two  cavities  are 
symmetrically  placed  within  the  cross  section  of  the  acceleration  unit,  the 
phase  shift  should  either  be  zero  or  it,  by  a  symmetry  argument.  The  case 
of  phase  shift  equal  to  n  is  equivalent  to  replacing  K  in  Eqs.  (2)  and  (3) 
by  -K.  Even  in  the  presence  of  the  beam,  we  have  found  from  the  numerical 
results  that  the  BBU  growth  is  not  sensitive  to  the  sign  of  K.  Henceforth, 
we  take  K  >  0. 

It  is  of  interest  to  note  at  this  point  that  the  presence  of  cross¬ 
coupling  leads  to  two  distinct  modes,  associated  with  the  ±  sign  in  Eq. 

(5).  Thus,  the  beam  needs  to  drive  two  modes  (instead  of  one  when  such  a 
coupling  is  absent).  For  the  same  beam  deflection,  the  deflecting  mode 
amplitude  in  a  cavity  would  be  less  when  k  ^  0.  This  reduction  in 
deflecting  mode  amplitude  would  in  turn  produce  less  transverse 


6 


displacement  the  beam.  This  is  believed  to  be  the  reason  why  cross¬ 
coupling  reduces  BBU  growth,  at  least  within  the  context  of  the  present 
model  of  cavity  coupling. 

Finally,  we  remark  that  the  "leakage"  of  mode  energy  to  another  cavity 
(when  K^)  is  not  equivalent  to  an  effective  lowering  of  Q  in  the  present 
model  of  cross-coupling.  A  finite  Q  always  represents  lossy  processes  from 
which  energy  cannot  be  recovered,  whereas  a  nonzero  K  represents  only 
reactive  loading  which  does  not  result  in  any  energy  loss.  This  is  obvious 
in  the  coupled  circuit  shown  in  Fig.  3.  To  see  this  mathematically,  we  set 
Q=®  (i.e.,  remove  the  lossy  processes)  but  keep  in  Eqs.  (2),  (3). 

These  two  equations  imply  conservation  of  total  mode  energy  because  of  the 
self-adjointness  of  the  operator  L  when  Q=®.  The  remaining  question,  which 
we  shall  not  address  in  this  paper,  is  whether  the  cross-coupling  is 
sufficiently  weak  to  allow  the  coupled  mode  treatment  given  in  this 
section.  The  alternative  would  be  to  treat  the  entire  accelerating  unit, 
including  all  of  the  gaps  on  the  various  arms  threading  through  the  unit, 
as  a  single  resonant  structure. 
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III.  Formulation  of  Beam  Breakup  Growth  in  the  Presence  of  Cavity 

Cross-Coupling 

The  formulation  of  beam-gap  interaction  follows  closely  those  given  in 
the  BBU  literature.  Consider  a  uniform  coasting  beam  of  pulse  length  T, 
relativistic  mass  factor  y.  and  current  I.  Let  x^^^(t)  and  p^^^(t),  be 
respectively,  the  transverse  displacement  and  the  transverse  momentum  of  a 
beam  slice  that,  at  time  t,  is  about  to  enter  the  i-th  accelerating  gap 
(i  =  1,  2,  ...20).  Let  x^^^(t)  and  p^^^(t)  be  the  corresponding  quantities 
when  this  beam  slice  exits  the  i-th  gap.  If  the  gap  width  is  small. 


^"^t)  =  x^^^t) 

(6) 

(i)(t)  =  p^^\t)  +  /^^t) 

(7) 

where  is  the  incremental  momentum  produced  by  the  deflecting  mode 

at  the  i-th  cavity  at  time  t.  The  i-th  cavity  is  in  turn  excited  by  the 
beam's  transverse  displacement,  and  because  of  the  cross-coupling,  it  is 
also  excited  by  the  neighboring  cavity  which  lies  on  the  same  plane 
orthogonal  to  the  beam  path.  Referring  to  Eqs.  (2)  and  (3),  the  evolution 
of  is  governed  by^’^ 

L/^'\t)  =  3  «^I  (^)  h^^^t)  x^^^t)  +  (8) 

c  o' 

i  =  1,  2,  ..20 

where  j  =  i+10  if  i  <  10  and  j  =  i-10  if  i  >  11. 

In  Eq.  (8),  (kc/w^)  is  a  dimensionless  quantity  equal  to  Zj^(C)/30  Q, 
and 


tj  <  t  <  t 

otherwise 


.  +T 
1 


(9) 
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represents  the  time  interval  during  which  the  beam  passes  through  the  1-th 
gap.  Ve  have  used  t^  to  designate  the  time  at  which  the  beam  head  enters 
the  i-th  gap.  If  T  designates  the  transit  time  between  neighboring 

f 

cavities  within  the  same  arm  and  T  the  transit  time  to  go  around  each  bend 
(e.g.,  from  cavity  no.  5  to  cavity  no.  6),  then  tj^  =  0,  ^2  ~  ^3  “  ^5 

=  4T,  tg  =  4T+t',  tj2  =  9T+2t',  etc.  [Fig.  2]. 

We  next  assume  that  the  beam  transport  between  neighboring  cavities 
within  the  same  arm  can  be  modeled  by  a  2x2  matrix  A.^  In  general,  a 

t 

different  matrix  A  would  be  necessary  to  model  beam  transport  across  the 
bend.^  Thus,  the  beam  transport  within  the  same  arm  is  described  by 


j^(i+l)^t+T) 

(t) 

=  A 

=  A 

pf  ^  (t). 

(10) 


for  i  ^  5,  10,  15.  Across  the  bends. 


x(l  +  l)(t4.T'  ) 

4-f 

x“>(,) 

=  A' 

II 

> 

^(i+l)^t+T'). 

(11) 


for  i  =  5,  10,  15.  In  writing  the  last  forms  of  Eqs.  (10)  and  (11),  we 
have  used  Eqs.  (6)  and  (7). 


Equations  (8),  (10)  and  (11)  are  the  governing  equations  for  the 
transverse  displacement  x^^^(t),  transverse  momentum  p^^\t),  and  the 
deflecting  mode  amplitude  /^^\t).  The  following  initial  conditions  and 
boundary  conditions  are  imposed.  We  assume  that  the  beam  enters  the  first 
cavity  without  any  initial  transverse  displacement  nor  transverse  momentum. 
For  the  cavities,  we  assume  that,  at  t=0,  they  are  all  unexcited,  except 
the  first  one,  which  has  an  amplitude  at  t=0.  Thus,  the  only  non¬ 
trivial  initial  condition  is 
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IV.  Nume^ic^l  Results 

In  the  presentation  of  the  numerical  results,  we  shall  use  the 
normalized  time  variables  t=(o^t.  Ve  shall  normalize  both  p  and  /  with 
respect  to  the  initial  amplitude  of  the  deflecting  mode  in  the  first 
cavity  [cf.  Eq.  (12),  Bq.  (7)].  The  diagonal  elements  of  the  transport 
matrices  A  and  A'  are  always  dimensionless.  The  off-diagonal  ones  become 
dimensionless  if  ve  normalize  p  by  Ve  shall  use  the  normalized 

variables,  unless  otherwise  specified.  The  normalized  equations  are  given 
in  Appendix  B. 

In  these  normalized  variables,  the  breakup  mode  period  in  the 

individual  cavity  is  2ii.  We  set  «  T  =  7.33,  «  T'  =  73.3,  »  t  =  204  in  all 

o  o  o 

numerical  calculations.  Note  that  the  pulse  length  is  sufficiently  short 
that,  at  any  instant,  only  one  (not  both)  of  the  cross-coupled  cavities  is 
occupied  by  the  beam.  The  above  numbers  correspond  to  a  breakup  mode 
frequency  of  1  GHz,  a  cavity  separation  within  each  arm  of  35  cm,  and  a 
distance  around  bend  (from  Gap  5  to  Gap  6)  of  350  cm.  To  characterize  the 
BBU  interaction  within  a  gap,  ve  introduce  a  dimensionless  parameter 


which  measures  the  strength  of  BBU  excitation  within  each  gap  [cf.  Eq.  (B3) 

of  Appendix  B].  For  I  =  10  kA,  y  =  17,  kc/w^  =  0.416,  ve  have  =  0.0144. 

This  value  of  e  has  been  used  in  all  numerical  calculations, 
r 

Finally,  for  the  transport  matrices,  we  assume  that,  for  simplicity, 

the  effective  betatron  frequencies  O)^  associated  with  beam  transport  within 

the  straight  arms  A,  B,  C,  D  are  identical  to  those  in  the  bends.  Ve  shall 

use  Q  ■  u  /w  to  denote  focal  strength.  Note  that  if  w  /2n  =  1  GHz,  0 
CO  o 

would  assume  a  value  of  0.28  if  the  effective  betatron  wavelength  2nc/u^  is 
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1.07  meter. 


L'or  Q  c  0.8,  the  elements  in  the  normalized  matrix  A  are 
®11  “  ®22  °  0'9134,  aj^2  “  -0.509,  32^  *  0.326,  whereas  those  in  A'  are: 
a'jj  =  a'22  =  -0.462,  a'j2  =  1.11*  2\  ~  -0.462.  (See  Appendix  BJ. 

In  the  numerical  calculations,  the  parameters  ve  have  used  are  0^4, 
20,  100,  10^;  S  =  0.14,  0.28,  0.8;  and  K  =  0,  0.01,  0.03,  0.1,  0.3.  These 
numbers  span  a  rather  vide  range.  [The  data  for  K  =  0.3  are  included  only 
to  show  the  qualitative  behavior  when  substantial  cavity  cross-coupling  is 
present.  One  may  argue  that  the  mode  coupling  formulation  is  adequate  only 
for  K  <  0.1] . 

In  order  to  read  the  figures  in  more  detail,  ve  make  a  few  comments 
about  the  dimensionless  time  coordinate  in  Figs.  4-9.  In  these  units,  the 
beam  enters  cavity  1  at  t=0  and  cavity  11  (the  cavity  to  which  it  is 
connected)  at  t*206.  It  enters  cavity  5  at  t=30,  cavity  6  at  t=103  (after 
a  bend),  cavity  16  at  t*310  (cavity  16  is  cross-coupled  to  cavity  6)  and 
cavity  20  at  t=340.  The  beam's  pulse  length  is  204  time  units. 

In  Figs.  4  and  5,  data  for  the  case  QslOO,  K=0  (no  cross-coupling), 
and  0=0.8  are  presented.  Figure  4  shows  the  deflecting  mode  amplitude  in 
cavities  1,  5,  6,  11,  16  and  20.  Figure  5  shows  the  beam  transverse 
displacement  in  the  same  cavities.  We  see  in  Fig.  4  that  the  deflecting 
mode  decays  continuously  with  time  in  cavity  1  as  expected  because  of  the 
finite  Q;  (cavity  1  is  not  excited  by  the  beam  since  the  beam  enters  on 
axis,  by  assumption.)  This  decay  continues  until  the  beam  has  left  the 
(un)coupled  cavity  11,  beyond  which  time  it  presents  no  physical  interest 
to  either  cavity  1  or  11.  The  field  in  cavity  5  shows  a  different  pattern. 
No  field  is  present  in  the  cavity  before  beam  entry  and  the  field  rises  in 
the  cavity  due  to  its  excitation  by  the  beam.  As  soon  as  the  beam  leaves 
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the  cavity,  tl^e  field  decays  as  in  cavity  1.  This  pattern  is  actually 
repeated  in  cavities  2  through  10.  The  pattern  in  cavities  11  through  20  is 
not  physically  different  but  appears  differently  in  the  figure.  In  these 
graphs,  only  the  rise  in  field  amplitude  is  shown,  since  the  decay  is  not 
relevant,  the  beam  not  re-entering  any  cavity  coupled  to  those.  Note  that 
in  this  case  of  no  cross-coupling,  the  field  amplitude  is  always  zero 
before  the  beam  enters  the  cavity,  except  in  cavity  1,  by  assumption. 

The  beam  transverse  displacement  does  not  show  as  much  structure  since 
the  beam  goes  only  once  through  each  cavity.  Figure  5  shows  the 
displacement  in  the  same  cavities  as  in  Fig.  4,  except  for  cavity  1,  which 
has  been  replaced  by  cavity  2  (the  displacement  in  cavity  1  is  always 
zero).  The  beam  picks  up  a  non-trivial  transverse  displacement  only  when 
it  arrives  at  cavity  2.  This  "initial"  displacement  is  equal  to 
x^=aj^2=®~^  sin(Qf)  [cf.  Eq.  (B-11)  of  Appendix  B].  In  order  to  assess  BBU 
growth,  we  take  the  ratio  of  the  final  beam  displacement  to  this  initial 
displacement  (Figs,  10,  11].  The  beam  displacement  in  Fig.  5  does  not 
evolve  exactly  as  the  field  amplitude  in  Fig.  4  but  it  can  be  seen  to  decay 
in  the  first  cavities  before  growing  exponentially  as  expected  of  BBU.  For 
the  present  case,  the  transverse  displacement  reaches  228  times  its  initial 
displacement  after  20  cavities. 

The  cavity  modes  shown  in  Fig.  6  correspond  to  the  same  Q,  and  Q  as 
in  Fig.  4  but  include  a  small  cross-coupling  K  =  0.01.  Although  the  final 
transverse  displacement  is  similar  (228  x^  vs.  210  x^),  the  mode  evolution 
is  qualitatively  very  different.  The  difference  starts  with  cavity  1,  for 
example,  where  the  field  decays  at  first.  However,  once  the  beam  enters 
the  coupled  cavity  11,  leakage  from  that  cavity,  although  not  very 
important,  shows  very  clearly  in  the  rise  of  the  field  in  cavity  1. 
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Although  this' field  will  not  act  on  2uiy  part  of  the  beam  directly  in  the 
present  configuration,  it  is  important  to  evaluate  it  since  it  affects  the 
field  in  cavity  11,  which  is  itself  driving  the  beam.  This  pattern  shows 
very  clearly  for  cavities  5  and  6.  In  cavity  11,  another  difference 
appears.  In  that  cavity,  when  the  beam  enters  it  at  tx206,  the  field  is 
non-zero  since  it  was  leaked  from  cavity  1  from  t=0  onward.  It  was  thought 
that  this  non-zero  field  could  aggravate  BBU  growth  but  it  turns  out  that 
the  opposite  conclusion  is  reached  -  namely,  that  this  non-zero  leakage 
actually  alleviates  BBU.  This  apparently  contradictory  result  will  be 
addressed  further  in  the  next  section. 

Differences  in  beam  behavior  are  much  smaller  than  differences  in 
field  amplitudes  and,  in  fact,  beam  displacement  is  quite  sirailar  to  the 
K=0  case  and  will  not  be  shown.  As  the  cross-coupling  constant  increases, 
the  coupling  between  modes  in  corresponding  cavities  increases  and,  as 
shown  by  relation  (5),  two  modes  would  be  excited  and  their  beating  give 
rise  to  much  more  complex  lime  behavior.  This  type  of  behavior  is  evident 
in  Fig.  7,  as  an  example,  where  K=0.1.  It  is  apparent,  by  comparing  Fig.  4 
and  Fig.  7,  that  increasing  the  coupling  coefficient  K  reduces  BBU  growth. 

Runs  have  been  repeated  for  the  case  Qs20,  0=0.28  [Fig.  8]  and  for 
0=4  and  0=0.14  [Fig.  9].  The  decreasing  values  of  0  correspond  to 
decreasing  values  of  the  focal  strength.  These  two  cases  display  behavior 
similar  to  the  higher  Q  case  shown  in  detail  previously.  The  maximum  beam 
transverse  displacement  at  cavity  20  normalized  to  its  initial  displacement 
is  shown  as  a  function  of  K  for  these  cases  in  Figs.  10  and  11.  The  effect 
of  cross-coupling  is  evident.  However,  it  becomes  less  important  as  Q 
decreases.  This  can  be  explained  easily  since  for  low-Q,  the  modes  are 
damped  much  more  effectively  and  that,  by  the  time  the  beam  enters  a 
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coupled  cavit;,  the  fields  in  that  cavity  vill  decay  to  a  low  level.  {See 
Fig.  9].  Basically,  a  low-Q  cavity  is  not  very  different  from  an  uncoupled 
one  and  this  fact  is  borne  out  by  the  results.  Figure  12  shows  the  peak 
mode  amplitudes  as  a  function  of  K,  in  several  cavities  for  the  case  Q^IOO. 

Ve  have  considered  three  different  values  of  Q  because,  although  the 
dominant  break-up  mode  may  have  a  lov-Q  value,  the  residual  deflecting 
modes  may  have  much  higher  Q  values  and  it  is  those  modes  that  are  most 
difficult  to  control.  However,  although  the  reduction  in  growth  is  larger 
for  higher-Q  cavities,  final  beam  displacement  is  still  larger  for  these 
modes.  There  is  substantial  reduction  in  the  transverse  displacement  by 
factors  of  ten  or  more,  when  the  coupling  constant  K  exceeds  a  few  per 
cent. 

Ve  summarize  the  main  features  in  the  above  analysis. 

1.  BBU  growth  decreases  with  increasing  cross-coupling  between 
cavities. 

2.  Significant  growth  reduction  occurs  when  k  >  a  few  per  cent. 

3.  The  reduction  in  BBU  amplitude  can  be  a  factor  of  10  or  more. 

4.  Substantial  growth  reduction  is  achieved  for  a  large  range  in 
K  values  and  so  should  be  readily  implemented  experimentally. 

These  features  are  rather  significant  because  they  show  a  possible 
reduction  in  BBU  growth  in  a  recirculation  system  which  is  not  available  in 
a  linac  geometry. 

Based  on  the  above  findings,  we  now  compare  BBU  growth  in  a  cyclic 

accelerator,  in  a  linac,  and  in  a  recirculating  linac  with  cross-coupling. 

Strictly  cyclic  systems  admit  unstable  normal  mode  solutions  (for  a  given 

9  10 

cavity)  with  simple  exponential  growth  in  time.  ’  Cumulative  BBU  in  a 

5  6  8 

linac  exponentiates  at  a  fractional  power  of  time  ’  ’  (<  1/2)  and  it 
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therefore  may»be  said  to  be  less  virulent.  Here,  we  have  shown  that  BBU 
growth  may  further  be  reduced  with  cross-coupling  of  cavities  in  a 
recirculating  geometry,  such  as  SLIA. 
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V.  Discussions 


The  numerical  results  based  on  the  rudimentary  theory  led  to  the 
following.  First,  the  proof-of-concept  experiments  in  SLIA  as  currently 
planned  are  not  likely  to  be  vulnerable  to  BBU.  Second,  contrary  to 
customary  thinking,  cross-coupling  of  cavities  in  different  arms  of  the 
SLIA  is  found  to  reduce  BBU  growth.  Since  the  second  point  is  rather 
unusual,  we  shall  have  to  explore  further  the  implications.  Me  shall  also 
comment  on  the  limitations  of  our  model  and  raise  the  various  issues  that 
need  to  be  addressed  if  SLIA  is  to  be  scaled  to  much  higher  energies. 

We  can  name  two  reasons  why  cross-coupling  among  cavities  could  reduce 
BBU  growth.  First,  when  coupling  is  present,  the  energy  in  one  cavity  is 
shared  with  another  cavity.  Thus,  the  cavity  does  not  provide  as  strong  a 
transverse  deflection  on  the  beam.  In  addition,  the  beam  needs  to  excite 
both  cavities  when  cross-coupling  is  present.  For  the  same  transverse 
displacement  of  the  beam,  the  deflecting  mode  that  would  be  excited  in  an 
individual  cavity  has  a  lower  amplitude.  Based  on  this  simple  physical 
argument,  it  is  natural  to  speculate  that  the  BBU  gain  per  pass  would  be 
further  reduced  if  each  accelerating  unit  were  to  accommodate  a  large 
number  of  arms.  In  such  a  case,  the  deflecting  mode  energy  in  one  cavity 
is  shared  by  a  large  number  of  other  cavities  through  mutual  cross¬ 
coupling.  This  scenario  is  quite  plausible  if  the  beam's  pulse  length  is 
sufficiently  short  that,  at  any  instant,  the  beam  resides  in  only  one  arm 
(among  the  various  arms)  in  the  cross-coupling  interaction.  Should  this  be 
the  case,  we  have  here  a  rather  peculiar  situation  for  SLIA;  The  BBU 
growth  per  pass  is  reduced  with  increasing  number  of  recirculations. 
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The  result  that  mode-coupling  can  reduce  BBU  growth  may  be  extended  to 
linacs.  If  the  accelerating  cavities  in  a  linac  are  coupled  to  similar, 
but  undriven  cavities,  BBU  growth  may  be  reduced  due  to  the  reactive 
loading  by  these  dummy  cavities.  [See  last  paragraph  in  Section  II]. 

Returning  to  the  present  study,  we  note  that  even  if  there  is  pre¬ 
existence  of  a  deflecting  mode  in  a  given  cavity  prior  to  beam  arrival  due 
to  cross-coupling,  this  initial  mode  energy  does  not  seem  to  reverse  the 
beneficial  effects  of  cross-coupling,  as  discussed  above.  Our  numerical 

4 

calculation  has  been  extended  to  very  high  Q  values  (Q=10  )  and  these 
conclusions  remain  valid.  Equation  (5)  suggests  the  existence  of  two 
distinct  modes  when  K^.  A  close  examination  of  the  numerical  results 
shows  that  the  solution  is  biased  toward  the  slow  mode  (o  =  a)^(l-K)  when  BBU 
is  excited. 

The  second,  but  perhaps  less  convincing,  reason  is  that  the  presence 

of  cavity  coupling  leads  to  a  detuning  of  the  breakup  mode  frequency  [cf. 

Eq.  (5)].  This  detuning  may  loosely  be  associated  with  an  intrinsic 

frequency  spread  which  is  known  to  provide  phase  mixing.  This  in  turn 

8  14 

leads  to  reduction  of  BBU  growth.  ' 

We  shall  now  address  the  deficiencies  in  our  model  and  speculate  on 
their  effects  on  the  BBU  growth.  Ve  have  assumed  that  the  values  of  Q  and 
associated  with  a  single,  isolated  cavity  remain  unchanged  when  this 
cavity  is  coupled  to  another.  The  adoption  of  this  assumption  may  have 
already  led  to  an  overestimate  of  the  BBU  growth  since,  intuitively,  when 
the  energy  within  a  cavity  is  allowed  to  couple  to  another,  Zj^  would  be 
reduced,  since  it  is  a  measure  of  the  excitability  of  the  cavity  by  a 
traveling  electron  that  is  displaced  from  the  center  axis  of  the  cavity. 
Other  simplifying  assumptions  ve  have  made  include  a  coasting  beam  and  the 
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use  of  constant  2x2  matrices  to  model  beam  transport  between  cavities 
(including  bends).  The  dependence  of  BBU  on  y»  which  varies  with  the 
focusing  system,  is  relegated  to  a  secondary  role. 

Let  us  now  consider  the  various  BBU  issues  that  need  to  be  addressed 
if  the  SLIA  is  extended  to  an  energy  much  beyond  8.5  MeV  (as  planned  in  the 
POCE).  In  that  case,  many  more  recirculations  and/or  higher  gain  per  arm 
would  be  required  and  each  beamlet  would  encounter  a  large  number  of 
cavities.  The  cross  section  of  an  acceleration  unit  would  accommodate  many 
arms  of  the  beam  pipe.  The  cross-couplings  among  these  various  arms  are 
clearly  far  more  complex,  though  they  may  lead  to  reduction  of  BBU  gain  per 
pass.  This  argument  actually  depends  on  whether  the  beam  consists  of  one 
single  pulse  or  a  series  of  pulses.  It  is  possible  that  our  conclusion 
regarding  BBU  reduction  by  cavity  cross-coupling  would  be  altered,  for 
example,  in  the  case  of  multiple  pulses.  One  can  think  of  situations  in 
which  cavity  cross-coupling  may  worsen  BBU  growth  in  a  recirculating 
geometry,  especially  when  the  individual  beam  pulses  simultaneously  occupy 
the  various  arms  within  the  same  acceleration  unit.  This  issue  remains  to 
be  investigated. 

As  we  have  emphasized,  SLIA  is  neither  a  linac  nor  a  strictly  cyclic 
system.  Applicable  scaling  laws  are  non-existent  at  the  time  of  this 
writing.  They  would  be  highly  desirable  but  their  derivation  touches  upon 
some  rather  complex  mathematical  issues.  Chief  among  them  is  the  question 
of  whether  conditions  exist  under  which  BBU  in  SLIA  would  evolve  in  a 
transient  manner,  or  into  an  unstable  normal  mode  (if  any). 

There  are  other  important  topics  that  have  not  been  addressed.  They 
are  mainly  related  to  ultra-high  energy  SLIA  operation.  For  example,  the 
numerous  encounters  of  the  beam  with  gaps  require  closer  attention  to  the 
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coupling  betveen  BBU  and  orbital  resonances . Perhaps  more  importantly 
the  many  recirculations  o£  a  high  current  beam  may  dangerously  couple  the 
negative  mass^^  and  image  displacement  instabilities  to  BBU.  These 
couplings  cannot  be  ignored  if  the  SLIA  is  to  be  extended  to  the  multi- 
kiloamp,  hundred  MeV  range. 
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Appendix  A.  Model  Equations  for  the  Coupled  Cavities 

Consider  first  the  coupling  of  the  two  identical  RLC  circuits  (labeled 
a,  b),  as  shown  in  Fig.  3.  These  two  circuits  are  linked  by  a  mutual 
inductance 

M  *  kL  (A1) 

where  K  <<  1  is  the  dimensional  coupling  constant.  Let  V  be  the  voltage 

across  the  inductance  L  of  circuit  a  and  be  the  current  flowing  through 

this  inductance.  Let  and  be  the  corresponding  quantities  for  the 

inductance  of  circuit  b.  Assume  solutions  with  exp(i(ot)  dependence.  The 


definition  of  mutual  inductance  M  gives 

V  =  i«LI  +i<oMI.  (A2) 

a  a  b 

=  i«LIjj+i«MI^  (A3) 

We  next  apply  the  Kirchoff  current  law  to  circuit  a  to  yield 

^  +  i«CV^  =  0.  (A4) 

Similarly,  for  circuit  b, 

V. 

Ijj+  ^  +  iwCVjj  =  0.  (A5) 


The  four  equations  (A2)-(A5)  in  four  unknowns  V^,  yield  the 

natural  frequencies  of  oscillation  of  the  coupled  circuits: 

(0  =  Wo[i+|q  ±  ^j.  (A6) 

2  2 

In  writing  (A6),  we  have  ignored  quantities  of  order  K  ard  1/Q  ,  and  we 

2 

have  used  the  relation  =  1/(LC)  and  Q  =  w^RC.  The  two  modes  given  in 
Eq.  (A6)  have  obvious  interpretation  in  the  Q  ->  ®  limit  and  in  the  K  ->  0 
limi t . 
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The  cross-coupling  between  the  cavities  shown  in  Fig.  2  may  now  be 
described,  making  use  of  the  preceding  results.  For  example,  let  circuit  a 
represent  cavity  no.  1  and  circuit  b  represent  cavity  no.  11  of  Fig.  3. 

The  field  amplitudes  and  may  then  be  described  by 


(A7) 

(A8) 

2  2  2 

where  L  s  d  /dt  +  <«^/Q)d/dt  +  Equations  (A7)  and  (A8)  are  adequate, 

as  far  as  mode  coupling  is  concerned,  because  the  two  normal  modes  given  by 
Eq.  (A6)  are  recovered  from  (A7)  and  (A8)  when  we  assume  normal  mode 

2 

solutions  of  the  form  exp(i(»>t),  under  the  same  conditions  K  «  1  and  1/Q 
«  1.  No  extraneous  modes  are  introduced,  and  Eqs.  (A7)  and  (A8)  also  have 
obvious  interpretations  in  the  limit  K  0. 
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Appendix  B.  Che  Normalized  Equations 

In  this  Appendix,  ve  normalize  the  governing  equations.  In  so  doing, 
we  identify  the  dimensionless  parameters  which  need  to  be  specified.  Ve 
use  a  bar  to  denote  normalized  quantities,  and  a  summary  is  given  at  the 
end  of  this  Appendix. 

Let  t  =  co^t,  f  =  w^T,  f'  =  «qT',  t  =  w^T,  f  s  f/f^, 
p  -  p//q»  X  =  yw^u^x/f^  where  is  the  initial  breakup  mode  amplitude  in 
the  first  cavity  [cf.  Eq.  (12)],  and  x,  T,  T' ,  are  respectively,  the  beam's 
pulse  length,  the  transit  time  between  neighboring  cavities  within  the  same 
arm,  and  the  transit  time  across  a  bend.  Ve  further  define 

L  B  d^/dt  ^  +  (l/Q)d/dt  +  1.  (Bl) 

Then  Eq.  (8)  becomes 

L  f  =  ej.h^^^t)x  ^^^(t)  +  K  7  (B2) 

where  i  =  1,  2,  3,  ...20;  j  *  i+10  if  i  <  10  and  j  =  i-10  if  i  >  11.  In 
Eq.  (B2),  h^^^  is  defined  in  Eq.  (9)  and  we  introduce  the  dimensionless 
constant  [cf.  Eq.  (13)] 


which  is  a  measure  of  the  BBU  strength  in  the  individual  beam-cavity 
interaction. 

The  beam  transport  Eq.  (10)  between  adjacent  cavities  of  the  same  arm 
is  then  normalized  to  read 


L(i^l)  _  _  1 

L  (i)  _ 

X  (t+T) 

A 

X  (t) 

s  A 

(i)  (i) 

Ip  (t)^/  (t)J 

i  /  5,  10,  15 


(B4) 
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where 


cos (ST)  Q~^sin(ST) 

-Qsin(OT)  cos(S7r)J 

and 

Q  E  u  /oD  (B6) 

c  o 

is  the  ratio  of  the  effective  betatron  frequency  (of  the  focusing  system 
within  the  arms)  to  the  breakup  mode  frequency.  Similarly,  the  transport 
equation  across  the  bend  [Eq.  (11)]  is  normalized: 


L  (i+l)  _ 

X  (t+T) 

-  A 

r  (i)  1 

X  (t) 

-(1+1)^-^^' 

s  A 

_(i)  _  _(i)  _ 

ip  (t)+/  (t)J 

i  =  5,  10,  15 

where 

cos(C'T)  (s')  sin(B'T) 
-Q'sin(Q'T  )  cos(fi'T  ) 

S'  E  w  '/(O 
c  o 


(B7) 


(B8) 


(B9) 


is  the  ratio  of  the  effective  betatron  frequency  (across  the  bend)  to  the 
breakup  mode  frequency. 

We  assume  initial  rest  conditions  for  all  and  except 

one,  namely, 

7^^^0)  =  1,  (BIO) 

as  in  Eq.  (12).  Thus,  the  beam  begins  to  pick  up  a  non-trivial  transverse 
displacement  only  when  it  arrives  at  cavity  no.  2.  Inserting  i  =  1  and  t  = 
0  in  Eq.  (B4),  we  obtain 

x^^^f)  =  aj27^^^0)  -  aj2  =  B"^sin(Bf)  e  (Bll) 
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which  is  the  normalized  tremsverse  displacement  of  the  beam  head  when  it 
first  encounters  the  second  cavity.  All  beam  transverse  displacement  in 
subsequent  times  should  be  compared  with  the  value  given  in  Eq.  (Bll)  to 
assess  the  BBU  growth. 

In  summary,  the  normalized  equations  for  cavity  excitation  are  given 
in  Eq.  (B2).  The  beam  transport  is  described  by  Eq.  (B4)  within  each  arm, 
and  by  Eq.  (B7)  across  the  bends.  The  only  non-trivial  initial  condition 
is  given  by  (BIO)  on  the  breakup  mode  amplitude  of  the  first  cavity. 
Besides  the  normalized  time  scales  x,  T,  T' ,  we  need  to  specify  e^,  0,  S', 
defined  in  Eqs.  (B3),  (B6),  (B9),  respectively.  Normalized  transverse 
displacement,  in  all  cavities  and  for  all  times,  should  be  compared  with 
the  value  given  in  Eq.  (Bll)  to  assess  the  increase  in  the  transverse 
displacement  as  a  result  of  BBU. 
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SLIA  PROOF-OF-CONCEPT  EXPERIMENT 


Accelerating  Units 


Fig.  1  —  Schematic  drawing  of  the  Spiral  Line  Induction  Accelerator 
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CROSS-COUPLING  OF  CAVITIES 


A 


D 


Fig.  2  —  Schematic  representation  of  cross-coupling  in  the  twenty  cavities 
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Fig.  4  —  Evolution  of  the  deflecting  modes  (i)  at  the  i-th  cavity.  In  the  left  column  (i  =  I,  S,  6).  the  arrows  on  the  left 
designate  the  time  at  which  the  entire  beam  exits  the  i-th  cavity;  the  arrows  on  the  nght  designate  the  time  when  the  whole 
beam  exits  the  (i-i-IO)th  cavity.  In  the  right  column  (i>=ll,  16,  20),  the  arrows  on  the  left  designate  the  time  of  arrival  of 
the  beam  head  at  the  i-th  cavity;  the  arrows  on  the  right  designate  the  time  when  the  entire  beam  clears  the  i-th  cavity.  The 
values  of  Q,  x,  and  Q  are  specified  on  the  figure.  The  coordinates  have  linear  scales. 
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Fig.  S  —  The  evolution  of  the  nonnalized  transverse  displacement  x  (i)  at  the  i-th  cavity.  The  arrows  in  the  left  column 
(i  =  2,  S.  6)  designate  the  time  at  which  the  entire  beam  leaves  the  i-th  cavity,  in  the  right  column  (i  =  Il.  16.  20).  the 
arrows  on  the  left  designate  the  time  of  arrival  of  the  beam  head  at  the  i-th  cavity;  the  arrows  on  the  right  designate  the  time 
at  which  (he  entire  beam  leaves  the  i-th  cavity. 
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Fig  6  —  Evolution  of  the  deflecting  modes  7  *  (r)  at  the  i-th  cavity.  In  the  left  column  (i=  1.  5,  6),  the  arrows  on  the  left 
designate  the  time  at  which  the  entire  beam  exits  the  i-th  cavity;  the  arrows  on  the  right  designate  the  time  when  the  whole 
beam  exits  the  (i-t-lO)th  cavity.  In  the  right  column  (i=ll.  16.  20),  the  arrows  on  the  left  designate  the  time  of  arrival  of 
the  beam  head  at  the  i-th  cavity;  the  arrows  on  the  right  designate  the  time  when  the  entire  beam  clears  the  i-th  cavity.  The 
values  of  Q,  v.  and  Q  are  specified  on  the  figure.  The  coordinates  have  linear  scales. 
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K  =0.1 

Q  =  100 

n  =  0.8 


Fig.  7  —  Evolution  of  the  deflecting  modes  /'*  (t)  at  the  i-ih  cavity.  In  the  left  column  (i  =  1,  5,  6),  the  arrows  on  the  left 
designate  the  time  at  which  the  entire  beam  exits  the  i-th  cavity;  the  arrows  on  the  right  designate  the  time  when  the  whole 
beam  exits  the  (i+IO)th  cavity.  In  the  right  column  (i=Il.  16.  20),  the  arrows  on  the  left  designate  the  time  of  arrival  of 
the  beam  head  at  the  i-th  cavity;  the  arrows  on  the  right  designate  the  time  when  the  entire  beam  clears  the  i-th  cavity.  The 
values  of  Q.  x,  and  Q  are  specified  on  the  figure.  The  Coordinates  have  linear  scales. 
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Fig.  8  —  Evolution  of  the  deflecting  modes  /  ’  (r)  at  the  i-th  cavity.  In  the  left  column  (i  =  I,  5,  6).  the  arrows  on  the  left 
designate  the  time  at  which  the  entire  beam  exits  the  i-th  cavity;  the  arrows  on  the  right  designate  the  time  when  the  whole 
beam  exits  the  (i-t-IO)th  cavity.  In  the  right  column  (i=ll,  16,  20),  the  arrows  on  the  left  designate  the  time  of  arrival  of 
the  beam  head  at  the  i-th  cavity;  the  arrows  on  the  right  designate  the  time  when  the  entire  beam  clears  the  i-th  cavity.  The 
values  of  Q.  »,  and  0  are  specified  on  the  figure.  The  coordinates  have  linear  scales. 
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Fig.  9  —  Evolution  of  the  deflecting  modes  (7)  at  the  i-th  cavity.  In  the  left  column  (i  =  1,  5.  6),  the  arrows  on  the  left 
designate  the  time  at  which  the  entire  beam  exits  the  i-th  cavity;  the  arrows  on  the  right  designate  the  time  when  the  whole 
beam  exits  the  {i-HO)th  cavity.  In  the  right  column  (i=Il,  16,  20).  the  arrows  on  the  left  designate  the  time  of  arrival  of 
the  beam  head  at  the  i-th  cavity;  the  arrows  on  the  right  designate  the  time  when  the  entire  beam  clears  the  i-th  cavity.  The 
values  of  Q.  x,  and  Q  are  specifled  on  the  figure.  The  coordinates  have  linear  scales. 
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Fig.  12  —  The  maximum  mode  amplitude  at  cavities  nos.  II,  16.  20  as  a  function  of  ».  Here  Q=  100,  0=0.8. 
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